Impurity-Induced Electronic Nematic State and 
C2-Symmetric Nanostructures in Iron-Pnictide Superconductors 
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We propose that impurity-induced electronic nematic state is realized above the orthorhombic 
structure transition temperature Ts in iron-pnictide superconductors. In the presence of strong 
orbital fluctuations near Ts, it is theoretically revealed that a single impurity induces local orbital 
order with C2-symmetry, consistently with recent STM/STS measurements. Each impurity-induced 
C2-symmetric nanostructure aligns along a-axis by applying tiny uniaxial pressure along 6-axis. In 
this impurity-induced nematic phase, the resistivity shows sizable in-plane anisotropy {pb/pa ~ 2) 
even above Ts, actually observed in various "detwinned" samples. The present study indicates the 
existence of strong orbital fluctuations in iron-pnictide superconductors. 

PACS numbers: 74.70.Xa, 74.20.-z, 74.20. Rp 
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I. INTRODUCTION 

Since the discovery of iron-pnictide superconductors,^, 
a lot of effort has been devoted to understand the overall 
phase diagram, including the superconducting (SC) state 
in the tetragonal (T) phase and non-SC orthorhombic 
(O) phase. In Ba(Fe,Co)2As2, the O structure transi- 
tion at Ts is second-order^, and very large softening of 
shear modulus Cs suggests the existence of strong ferro- 
quadrupole 0s oc — fluctuations above Ts'^"^. In 
many compounds, the superconducting (SC) transition 
temperature Tc takes the highest value near the endpoint 
of the O phase, suggesting a close relation between the 
superconductivity and the orbital instability. The weak 
ferro-quadrupole order in the O phase induces the spin- 
density- wave with Q = {it, 0)^. 

As for the SC mechanism, spin fluctuation mediated 
sign-reversing s-wave state (s±-wave state) had been pro- 
posed from the early stage, noticing on the Coulomb 
interaction and intra-orbital nesting between hole- and 
electron-pockets^^^. In iron-pnictides, each Fermi pock- 
ets are mainly composed of t2g orbitals of Fe atoms. 
On the other hand, orbital fluctuation mediated s-wave 
state without sign reversal (s+^-wave state) had been 
investigated in Refs.^""^^: Strong orbital fluctuations 
originate from the inter-orbital nesting between Fermi 
pockets in the presence of Coulomb and weak electron- 
phonon (e-ph) interactions. The latter scenario is sup- 
ported by the robustness of the SC state against im- 
purities in many iron-pnictides"'^^^^^, and by the or- 
bital independent SC gap observed in Bal22 systems 
by laser ARPES measurement^^'^^ Also, experimental 
"resonance-like" hump structure in the neutron inelastic 
scattering is well reproduced in terms of the s++-wave 
SC state, rather than the s±-wave SC state, by taking 
the suppression in the inelastic scattering in the SC state 
(dissipationless mechanism)^*. 

According to Ref.^, the structure transition originates 
from the ferro-charge quadrupole (j)s — 0.^.2 _y2 oc rixz — 
Uyz instability, realized by the bound state formation 
of two orbitons with opposite momenta. By this two- 



orbiton theory, we can fit the temperature dependence 
of Cs in Ba(Fei_a;Coa;)2As2 for a; = ~ 0.16 almost 
perfectly The spin nematic theory (or two-magnon 
process)^'^" is another candidate. However, incommen- 
surate spin order is realized in Ba(Fei_a;Coa;)2As2 for 
X > 0.056^^, although the latter theory requires com- 
mensurate fluctuations. 

Furthermore, recent discovery of "electronic nematic 
transition" in the T phase, free from any lattice deforma- 
tion, has been attracting great attention. For example, 
in "detwinned" Ba(Fei_a;Coa;)2As2^^'^'^ under very small 
uniaxial pressure (~ 5MPa), sizable in-plane anisotropy 
of resistivity emerges at T*, which is about lOK^lOOK 
higher than Ts- The nematic order is also observed in 
BaFe2(As,P)2 by the magnetic torque measurement^*. 
Now, it is demanded to flnd which degree of freedom, 
spin or orbital, is more important for the nematicity, or- 
thorhombicity and superconductivity. 

In this paper, we discuss the impurity-induced elec- 
tronic nematic phase in iron-pnictides, using the mean- 
fleld approximation (MFA) in real space. When or- 
bital fluctuations develop, we obtain various types of 
local orbital orders with lower symmetries (C4, C'2v, 
C2, etc.), actually reported by STM/STS autocorrela- 
tion analyses^^'^^. The large cross section of the local 
order gives giant residual resistivity, far beyond the s- 
wave unitary scattering value; 20/if2cm/%. When 
C2 nanostructures are aligned along a-axis, the in-plane 
anisotropy of resistivity reaches 40%, consistently with 
cxperiments^^'^^. Such large anisotropy is not achieved 
when isotropic impurity scattering is considered^ 

In annealed Ba(Fei_^M2;)2As2 (M=Co, Ni), the differ- 
ence Ipfe — Pal is very small in the absence of M-impurities 
{x = 0)^^, while it increases in proportional to x for 



< 4% 



23,29 



In contrast, both the magnetic moment 



and lattice deformation monotonically decrease with x. 
These facts strongly support the idea of impurity-induced 
nanostructures . 

In strongly correlated electron systems, impurity po- 
tential frequently causes drastic change in the electronic 
state. For example, in nearly antiferromagnetic metals. 
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magnetic correlation is extremely enhanced near the non- 
magnetic impurity site, giving rise to the local magnetic 
moment I/xb) and large residual resistivity ^° that are 
indeed observed in optimally- and under-doped cuprates. 
In terms of weak-coupling scheme, such phenomena orig- 
inate from the Friedel oscillation since the large local- 
density-of-states (LDOS) sites could trigger the strong 
fluctuations around the impurity. As for the iron-based 
superconductors, the system would be close to antiferro- 
orbital critical point. Thus, it is natural to expect the 
occurrence of "impurity-induced local orbital order" in 
iron pnictides. 



II. MODEL HAMILTONIAN AND METHOD OF 
CALCULATION 

Here, we study the single-impurity problem due to 
orbital-diagonal impurity potential I^" in a large clus- 
ter with 800 Fe sites, based on the two-dimensional ten- 
orbital tight-binding model for LaFeAsO in Refs.^'^^. We 
set X and y axes parallel to the nearest Fe-Fe bonds. 
Then, the Fermi surfaces are mainly composed of t2g or- 
bitals {xz, yz and xy), although Cg orbitals also play 
non-negligible roles. Here, we consider the following 
quadrupole-quadrupole interaction^^^""^^: 



quad 



XZ XZ 



o: 



6' \ (V) 



where Op is the quadrupole operator for channel F at site 
i introduced in Ref.'^: Of = J2i,m,a Or'^^l,icr'^i,ma, where 
Op™ is the matrix element of the charge quadrupole op- 
erator. Note that 0^^ oc lfj_li, + l^l^- The quadrupole 
coupling constant g in eq. (1) originates from both the 
e-ph interaction as well as the Coulomb interaction for 
the charge sector, as discussed in Refs.^"'^^. Since we 
are interested in the nonmagnetic orbital order, we ne- 
glect the Coulomb interaction to simplify the calculation. 
Then, strong orbital fluctuations for F = xz,yz channels 
are produced by relatively small g (~ 0.2 eV)^*^"^^. Here- 
after, the unit of energy is eV. 

Here, we put T = 0.02 and the electron filling n = 6.0 
per Fe, which corresponds to undoped compounds like 
BaFe2As2. In the absence of impurity, the bulk antiferro- 
orbital order occurs for g > gc = 0.222. Below, we study 
the following mean- field equation for g < g^- 



Ml 



/,0 
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where i is the Fe site, and l,m represent the d-orbital. 
M/^ is impurity- induced mean-field; = for / = 0. 
Then, the mean-field potential due to Hartree term is 



hn.V 
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where Ff^ 



xz,yz,xy Am J' m' 



r 



is the bare inter- 



is Hmf = Hq + Y^-S'' + const. 
(2)- (3) self-consistently. 



In the MFA, we solve eqs. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 
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action for charge sector^, and the mean-field Hamiltonian 



FIG. 1: (color online) Obtained LDOS at Ef for I = -2 and 

(a) g = 0.200: without orbital order, (b) g = 0.208: orbital 
order with diagonal C2v -symmetry, and (c) g = 0.218: orbital 
with C2-symmetry. (d) Energy-dependence of the LDOS for 
g — 0.218. (e) g-dependence of the free-energy. 

In Figs. 1 (a)-(c), we show the obtained DOS at Fermi 
level (E-p) in real space, in which the center is the impu- 
rity site with / = —2. For g = 0.200 (a), the impurity- 
induced mean-field is absent. The small modulation of 
the LDOS around the impurity is caused by the Friedel 
oscillation. For g > 0.207, impurity- induced local orbital 
order with diagonal C2v symmetry appears, as shown in 

(b) . The suppression of the DOS is caused by the or- 
bital order, consistently with a recent optical conductiv- 
ity measurement^^. With increasing g, the orbital order 
changes to C2 symmetry for g > 0.212, shown in (c). 
The size of the nanostructure is ^ 15aFc-Fc ^apc-Fc) 
along x (y) axis. Such a large impurity-induced object 
is actually observed in Ba(Fe,Co)2As2 by STM/STS25'26. 

(3) When the impurity concentration njmp is 1%, the ob- 
tained C2-order would be stabilized by the weak overlap 
between neighbors against thermal fluctuations omitted 
in the MFA. (Similar C2-order is also realized for / = 00.) 
Figure 1 (d) shows the energy-dependence of LDOS for 
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g = 0.218 at r = (0,0) (impurity site), (1,0), (2,0), 
and (4,4). Near (0,0), the LDOS is modified for a wide 
energy range. Figure 1 (e) presents the free-energy as 
function of g. In the MFA, each transition dX g k, 0.207 
and 0.212 is first-order. 



n 



100 




-6 -4 -2 2 4 6 -6 -4 -2 2 




-6 -4 -2 2 4 6 -6 -4 -2 2 4 6 



FIG. 2: (color online) Obtained electron density n' and 
quadrupole order Or at Fe sites for I = —2 and g = 0.218: 



Since M, 



l,m 



J, the present mean- field has 15 com- 



ponents at each site. They are represented as charge den- 
sity or monopole (1 = 0), quadrupole (/ = 2), and hex- 
adecapole {I — 4) orders. The first two orders are give 
as = 2Em^^m and Of = 2 E^,™ o[;™M/,„, where 
r — xz, yz, xy, z'^, and x'^ — y^. The hexadecapole order 
is negligibly small in the present study. Figure 2 shows 
the dominant four mean- fields, n*, O^^, Oy^ and O^y, 
for g = 0.218. We verified that the quadrupole interac- 
tions for F = xz/yz channels in eq. (1) are indispens- 
able for the C2-order. The obtained quadrupole order 
is very difference from the uniform quadrupole ordered 
state (0^2_y2 oc Uxz — n-yz =const.) in the orthorhom- 
bic phase^, and therefore the impurity-induced nematic 
order will exist even below T5. 

The C2-order in Fig. 1 (c) can be aligned by the strain- 
induced quadrupole potential; H' — AE ^ - OI2 .,2 and 



(0)/xil ,2(0), where eg 



oc a 



AE = rjses ■ X^2„y. 
is the strain and rjs is the strain-quadrupole coupling. 
Xx'^-y2{0) is the ferro-quadrupole susceptibility, which is 
strongly enhanced near T5 due to the two-orbiton process 
as discussed in Ref.^. This would be the reason why 
the nematic ordered state is easily detwinned by small 
uniaxial pressure near Ts- In fact, de twinning by uniaxial 
pressure is possible only when the structure transition is 
the second-order^'^. 

Here, we assume x \\ a-axis and y \\ 6-axis. In de- 
twinned compounds with a > b, ARPES measurements 
indicates AE < 0, i.e., Uxz > nyz^"^'^^ . For a single C2- 
order, we obtain the relation Fa — Ft, k, 2.5AE, where 
-Pa(6) is the free-energy when the C2-order is along a{b)- 
axis. Therefore, the nematic order along a-axis is realized 
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FIG. 3: (color online) (a) Alignment of the impurity-induced 
C2-orders under uniaxial pressure (a > 6). (b) Obtained pa(b) 
for riimp ~ 1% and / = —2: pb > pa in the nematic phase. 



by detwinning (a > b), schematically shown in Fig. 3 (a). 
Note that two kinds of C2-orders, the C2-order in Fig 1 
(c) and its inversion with respect to x axis, still degener- 
ate and coexist with equal probability. 

Now, we calculate the in-plane resistivity in the ne- 
matic state shown in 3 (a). We use the T- matrix approx- 
imation, which gives the exact result when riimp <C 1 and 
localization is negligible. The T-matrix is given by solv- 
ing the following equation in the orbital-diagonal basis: 

Tr,r'i^) = {I + S)rSr,r' 

+ Y,ii + S)rG^X"i^)Tr",r'ii0) (4) 
r" 

where Sr is the impurity-induced mean-field potential, 
and G^"^ (lu) is the Green function without impurities. 
Ir ~ liSr.o is the impurity potential^*'. The T-matrix 
is non-local when Sr ^ 0. After the Fourier transfor- 
mation, the self-energy in the T-matrix approximation 
is E(fc,a;) = ?iimpT'fc,fc(w), and the full Green function is 
G{k,uj) = {w + n - H° - T,{k,uj))-^ . Then, the in-plane 
conductivity is given as 



1 ^ 



(5) 



where v — x or y, and a represents the ath band, ^ 
is the group velocity and Ga{k,uj) is the full Green func- 
tion in the band-diagonal basis. ^ is the total current 
including the vertex correction, which is given by solv- 
ing the following Bethe-Salpeter equation: ^ = t'fe ,y + 

^E4%^|G^(P'*^)l'^P^^here4";f,=ni„,p|T,";,^,(z<5)p 

is the irreducible vertex. 

The obtained results for I — —2 and riimp = 1% are 
shown in Fig. 3 (b). Here, we assume the inter-layer 
distance is 0.6nm. Without orbital order, the resistiv- 
ity is 5.5/ir2cm, which is about one- fourth of the max- 
imum value without orbital order: puni 20^QcTa for 
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/ w +1. When diagonal C2i,-order appears, the resistiv- 
ity exceeds puni, due to large cross section of the "effec- 
tive impurity radius" as recognized in Fig. 1 (b). In the 
nematic phase with horizontal C2-order, we obtain large 
anisotropy Pb/ Pa ~ 2: By including the vertex correc- 
tion, both Pa and pb are suppressed and the anisotropy 
Pb/ Pa is enlarged, since the contribution of the forward 
scattering is correctly subtracted. The averaged resis- 
tivity {pa + Pb)!'^ per 1% impurity reaches ^ 50^f7cm, 
which is comparable to the residual resistivity by 1% Co 
impurities observed in Lallll^^ and Bal22^^. 

Now, we discuss the nematic transition at T* in real 
compounds. Beyond the MFA, the effective interac- 
tion g (< g) decreases with T due to the thermal 
fluctuation^^. Then, one possibility is that the phase 
transition from the diagonal to vertical Ci occurs 
at T*. (Then, ~g w 0.212 at T*.) Other possibility is 
that Ci order is realized even above T* , while the neces- 
sary condition for detwinning, x^2_j,2(0) ^ 1, is satisfied 
only below T*. In both cases, experimental results can 
be explained. 




I ■ U.Z U.Z I 

-4 -3-2 -1012 3 4 g [eV] 



FIG. 4: (color online) Obtained LDOS at .Ef in the case of 
/ = -1-1 and (a) g = 0.200: without orbital order, and (b) 
g = 0.218: orbital order with C4 symmetry, (c) Obtained 
resistivity for 1% impurity with / = +1. 

We also study the impurity-induced local orbital order 
for / — -1-1. Figure 4 (a) shows the LDOS without orbital 
order: The realized Friedel oscillation pattern different 
from Fig. 1 (a) would induce a new type of orbital order. 
In fact, we obtain the orbital order with C4 symmetry for 
g > 0.203: Figure 4 (b) shows the LDOS for g = 0.218. 
We also obtain a meta-stable solution with C2 symmetry 
similar to Fig. 1 (c), whose free-energy is about 0.1 eV 



higher than that for the C4 symmetry solution. (When 
/ = —2, the C4 symmetry solution is "unstable" with 
positive free-energy.) Figure 4 (c) shows the resistivity 
p — Pa — Pb ioT I = +1 and nimp — 1%: It exceeds the 
unitary value as soon as C4-order appears, and it reaches 
~ SO/ifJcm for 5 ge- 
lt is noteworthy that the obtained C4 order looks 
similar to Sn-impurity-induced "ring-shape object" in 
LiFeAs observed by Hanaguri"^^ very recently. The 
realized large reduction in the DOS would result in 
the suppression of the s.4__|_-wave state. In fact, in 
BaFei.89-2a:Zn2a;Coo.iiAs2, the suppression in Tc per 1% 
Zn-impurity is —ATc/% ^ 3 K/%^^: Such smaU suppres- 
sion of Tc is consistent with the s++-wave state, since 
— ATc/% ^ 20 K/% is expected in the s±-wave state 
when the mass-enhancement is m* /irib ~ 3^^. 

Finally, we consider the impurities on other than Fe 
sites. We expect that impurity-induced nematic phase 
is realized in BaFe2(As,P)2, since P sites give finite im- 
purity potential on the neighboring four Fe sites. In this 
case, we actually obtain impurity-induced order with C2- 
or Ci/i-symmetry, consistently with experiments^^. In 
contrast, nematic state is not realized in (K,Ba)Fe2As2'^^, 
maybe because K sites are outside of FeAs planes. 



IV. SUMMARY 

In summary, we discussed impurity-induced electronic 
nematic state based on the orbital fluctuation theory. 
The obtained local orbital orders with various sym- 
metries (C2V, C'2, and C4) are consistent with recent 
STM/STS measurements. In the case of C2-ordcr, 
the anisotropy of resistivity reaches Pb/pa ~ 2, which 
presents a natural explanation for the nematic state in 
various "detwinned" iron-pnictides. Thus, characteristic 
features of iron pnictides, nematic and structure transi- 
tions as well as superconductivity, are well understood 
based on the orbital fluctuation theory. 
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